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It has recently been proposed that Na4lr308 is a weak Mott insulator at ambient pressure, sup- 
porting a three-dimensional spin liquid phase with a spinon Fermi surface. This proposal is consistent 
with recent experimental findings that the material becomes a metal upon increasing pressure or 
doping. In this work, we investigate the effect of the spin-orbit coupling arising from 5d Ir moments 
both in the metallic and spin liquid phases of Na4lr308. The effective Hubbard model in terms 
of pseudospin j — 1/2 Ir states is derived and its consequences to both metallic and spin liquid 
phases are studied. In particular, the model leads to enhanced Wilson ratio and strong temperature 
dependence of the Hall coefScient. 



I. INTRODUCTION 

Recently, Na4lr308 has emerged as a promising candi- 
date for a three dimensional quantum spin-liquid. Sus- 
ceptibility measurements show a large antiferromagnetic 
Curie- Weiss temperature in this material, Qcw = —650 
K, yet no evidence of magnetic ordering has been ob- 
served down to very low temperatures^. The material at 
ambient pressure is an insulator, yet specific heat mea- 
surements support the existence of a finite density of low 
energy excitations. This has lead to different theoret- 
ical proposals for Mott insulator states with fermionic 
spinon excitations^"—, and to proposals to detect the 
spinon Fermi surface^. 

At a superficial level, a three dimensional material 
may be an unlikely candidate for a quantum spin-liquid. 
Fluctuations are usually weak in three dimensions, and 
magnetic order is likely to set in at finite temperatures. 
However, despite its three-dimensional nature, Na4lr308 
combines a confluence of factors which favour the for- 
mation of a spin liquid: the Ir atoms have spin-1/2, 
leading to strong quantum fluctuations; these atoms 
are arranged in an undistortcd "hyperkagome" lattice, 
which gives rise to magnetic frustration; and, perhaps 
most importantly, there are strong charge fluctuations 
in Na4lr308. These charge fluctuations lead to ring- 
exchangc-type interactions of increasingly long-range be- 
tween moments, and can prevent the formation of long- 
range order altogether—. 

The strongest evidence for charge fluctuations in 
Na4lr308 comes from recent measurements on samples 
under pressure and on doped samples^^. Application of 
hydrostatic pressure enhances the conductivity of the un- 
doped samples by several orders of magnitude, indicating 
the proximity of Na4lr308 to a metallic phase. This is 
further corroborated by the fact that relatively small con- 
centrations of dopants turn the system metallioi^. This 
indicates that Na4lr308 is a weak Mott insulator, which 
can readily undergo bandwidth and doping-controlled 
transitions to the nearby metallic phase. 

The weak Mott insulator scenario for Na4lr308 was 
first proposed in Ref. 0, in which the bandwidth-tuned 



transition to a metal was studied in detail. In this the- 
ory, the electrons are factorized into fermionic spinous, 
which carry spin but no charge, and bosonic rotors, which 
carry charge but no spiniS. When a critical interaction 
strength is reached, the rotors become gapped, and the 
system becomes an insulator. On the other hand, the 
spinous remain gapless even in the insulator, forming a 
neutral spinon Fermi sea, responsible for the finite den- 
sity of excitations seen in the specific heat. The spin- 
liquid states proposed so far, however, have ignored the 
sizeable spin-orbit coupling effects present in Na4lr303^. 
The heavy Ir atoms have a large spin-orbit coupling 
which must be included in the same footing as electronic 
hopping and correlation effects in any microscopic theory 
of Na4lr308. 

In this paper, we study the effects of spin-orbit cou- 
pling on the weak Mott insulator and metallic phases of 
Na4lr308. We find that the system is well-described by a 
one-band Hubbard model, in which the electrons occupy 
total angular momentum j ~ 1/2 ("pseudospin") states 
on sites of the hyperkagome lattice. Thus, the situation 
is similar to the absence of spin-orbit coupling, but with 
two new ingredients: instead of carrying spin s ~ 1/2, 
the electrons carry pseudospin j = 1/2, and the hopping 
between sites does not necessarily conserve pseudospin. 
However, the fact that the system is still described by a 
one-band model is significant, as it justifies the general 
approach followed in Ref. 0. 

The importance of spin-orbit in Na4lr308 can be ap- 
preciated from recent studies on other transition metal 
oxides containing iridiuroii"— . For example, one may ex- 
pect that correlation effects in Sr„+iIr„03„-)-i would be 
weak due to the extended nature of the 5d Ir orbitals. In- 
stead, spin orbit coupling gives rise to a narrow j ~ 1/2 
band in these materialsi^. In Sr2lr04 this leads to an 
antiferromagnetic Mott state, as demonstrated through 
angle resolved photocmission spectroscopy, optical con- 
ductivity, and x-ray absorption spectroscopy measure- 
ments, and also through first-principles electronic struc- 
ture calculationsi^ii^. We find that a similar mecha- 
nism is responsible for enhancing correlation effects in 
Na4lr308. 
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The inclusion of spin-orbit coupling may also help ex- 
plain one of the outstanding experimental observations 
in Na4lr308- At low temperatures in the Mott insulator 
phase, the specific heat has a linear temperature depen- 
dence, with a finite coefficient 7 = C/T, as in a metal. 
Similarly, the magnetic susceptibility is constant at low 
temperatures, as in a metal. These results may be evi- 
dence for the existence of a spinon Fermi surface in the 
insulator—. On the other hand, in the experiments the 
Wilson ratio between the susceptibility and the specific 
heat coefficient is much greater than one. This is of- 
ten taken to be a signature of strong ferromagnetic fluc- 
tuations, but it seems to be at odds with the antifer- 
romagnetic Curie- Weiss temperature^. The inclusion of 
spin-orbit coupling gives an alternate way of generating 
a Wilson ratio different from one^, which we will explore 
below. 

Our approach in this paper is partially phcnomcnolog- 
ical. After deriving an effective Hubbard model based on 
microscopic considerations, we are left with a single free 
parameter 6 that controls the ratio of direct Ir-Ir hopping 
to the hopping mediated by oxygens. Since the outer 5d 
electrons in Ir are extended and it is difficult a priori to 
make an estimate of 0, instead we follow a phenomenolog- 
ical approach and estimate 9 by requiring a large Wilson 
ratio and a small 7. However, many of our qualitative re- 
sults are independent of the precise value 0. For instance, 
we find that for all values of 9, the non- interacting system 
has a finite density of states. This is in contrast to the 
pyrochlore lattice, where spin-orbit effects can lead to the 
formation of a topological band insulator—. Therefore, 
the insulating behaviour in Na4lr308 must arise due to 
interaction effects. 

Another generic feature that is independent of the de- 
tailed value of is a strong temperature dependence in 
many physical quantities. The hyperkagome lattice has a 
12 site unit cell which, together with the presence of spin- 
orbit coupling, leads to 24 quasiparticle bands. We find 
that the Fermi surface contains multiple pockets arising 
from many of the bands, which have very strong struc- 
ture near the chemical potential. This leads to a strong 
dependence of the density of states with energy, which 
is reflected in the temperature dependence, 6.17. of the 
speciflc heat and also of the Hall coefficient for doped 
samples. 

This paper is organized as follows. In Section we 
give a microscopic derivation of the one band Hubbard 
model of Na4lr308in terms of the pseudospin j = 1/2 Ir 
states. In Section HI, we analyze the quasiparticle spec- 
trum in the non-interacting limit and study the resulting 
thermodynamics and transport properties in the metal- 
lic phase. We discuss the implications of the spin-orbit 
coupling in the spin liquid insulator. 



II. HUBBARD MODEL 

Our starting point in deriving a microscopic model for 
Na4lr308 is the observation that due to the large atomic 
number of Ir, the spin-orbit coupling Ago at the Ir sites is 
expected to be comparable to other microscopic energy 
scales, such as the local Hubbard repulsion U and the 
hopping amplitude between neighbouring Ir sites, and 
therefore must be treated on the same footing. Therefore, 
electronic spin 5' is not expected to be a good quantum 
number. In what follows, we first consider a single Ir 
site with spin-orbit coupling, which wc will then use to 
construct a tight binding model. 

The Ir 5d orbitals arc split into Cg and t2g levels by 
the octahedral crystal field of the nearby oxygens. The 
Sg — t2g splitting is large relative to A^d^ '^^'^"^ . Therefore, 
the Cg orbitals in Ir*"*" are completely empty and can be 
ignored, whereas the t2g orbitals have a single s = 1/2 
hole. The t2g orbitals are weakly split due to weak dis- 
tortions of the oxygen octahedra. However, wc assume 
that this splitting is much smaller than Ago- Then, the 
three t2g orbitals behave as an effective orbital angular 
momentum L ~ 1 moment. The spin-orbit coupling be- 
tween this L ~ \ moment and the s = 1/2 spin yields a 
low energy j = 3/2 multiplet and a high energy j = 1/2 
multiplct, separated by the energy Ago- Provided that 
the electronic hopping between Ir atoms is not large rel- 
ative to Aso, the j = 3/2 and j = 1/2 bands remain 
well-separated. Then, the j ~ 3/2 band is completely 
full and the j = 1/2 band is the only one to participate 
in the low energy physics. The emergence of a j = 1/2 
band in an iridium oxide has been studied in detail in the 
closely related material Sr2lr0 4-^^i-^^i-^^ 

Hence, j = 1/2 states give a good description of Ir elec- 
trons in Na4lr308. The j = 1/2 states are represented in 
the jz = ±1/2 basis by?^ 

where we have used the subscripts j and s to distinguish 
between total and spin angular momenta labels. By com- 
parison, in the absence of spin-orbit coupling, the de- 
scription would be in terms of s = 1/2 moments. In 
what follows, we need to keep track of whether the la- 
bels I t) a-nd I i) represent spin angular momentum s, 
or "pseudospin" j. We will use the labels a = ±1 and 
a ~ ±1 to describe pseudospin and spin, respectively, i.e. 

I« = +1) = 1 1,), l« = -1) - I i,), k = +1) = I ts), 

and |a = — 1) = I ^s). We can now introduce interac- 
tions and hopping, to obtain a one-band Hubbard model 
on the hyperkagome lattice: 

i i 

Here, dia annihilates a d electron on the Ir site z, rii = 
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FIG. 1; (Color online) Neighbourhood of one Ir site (C) in 
the hyperkagome compound Na4lr308. Site C is surrounded 
by an octahedron of oxygens (labelled 1 ... 6), and it has four 
nearest neighbour Ir sites, A, B, D, and E. We ignore the 
weak cubic distortions present in Na4lr308. Note that site 
A has only one second neighbour shown in the figure (D) 
and one third neighbour (E), whereas site B has two second 
neighbours shown (D and E), and no third neighbours. 

J2a d'iadia IS thc total dcnslty of electrons occupying 
the j = 1/2 multiplet on that site, and Hq is a non- 
interacting tight-binding Haniiltonian, 

i?0 = - ^ tia.ja'dl^dja' ■ (2) 
ijaa' 

The pseudospin-dependent hopping amplitude tia^ia' is 
derived below. 



A. Tight-binding parameters 

Once we consider many sites, the hopping ampli- 
tudes become pseudospin-dependent due to spin-orbit 
coupling. To derive the effective tight-binding Hamil- 
tonian for the j = 1/2 bands, we assume three types 
of hopping processes are present: (i) Direct hopping be- 
tween d-orbitals on the nearest neighbour Ir sites {e.g. 
sites A and C in Fig. [TJ 

+t2d\^y^^„dc,xz,a + h.c. (3) 

(ii) Hopping between Ir c?-orbitals and the nearest neigh- 
bour O p-orbitals {e.g. between sites A and 1): 

HaI = -Upd\.^y^„Plxa - Upd\^y^^„Plza + h.c. (4) 

and (in) Hopping between nearest neighbour O p-orbitals 
{e.g. between sites 1 and 5): 

Hiz = -Up[^^p5,^ [p[y^P5ya + P[xaP5xa ) 

- 2 {p\yaP5xa + p\y^P5xc7^ + h.C. (5) 



In these expressions, summation over repeated spin in- 
dices a is implicit. Note that we have defined 5 indepen- 
dent hopping parameters, ti, t2, tdp, t^, and to-, where 
the last two describe tt and cr hopping processes between 
nearest-neighbour p-orbitals . 

Most of the hopping processes relevant to Na4lr308 
involve virtual states with vacancies in the oxygen p- 
orbitals. These processes are suppressed by the energy 
denominator Ep = ed — — Up > 0, where and ep are 
the level energies of the j = 1/2 iridium orbitals and the 
oxygen p orbitals, respectively, and Up is the Hubbard re- 
pulsion on the oxygens. The energy Ep is assumed to be 
large relative to the hopping amplitude between iridium 
and oxygen, and also to the hopping between two oxy- 
gens. Therefore, the energy Ep controls the perturbation 
theory used to compute the effective tight binding model 
parameters, which are expressed in powers of t/Ep. 

Symmetry arguments provide strong constraints on the 
form of the tight-binding model that obtains from per- 
turbation theory. For instance, since spin-orbit coupling 
is time-reversal invariant, the most general pseudospin- 
dependent hopping appearing in Eq. ([2|) is of the form 

where v^j is a real vector and a = {<Jx,'^y,'^z) arc the 
three pseudospin Pauli a matrices. Note that Vy must 
be odd under exchange of i and j in order to ensure 
Hermiticity of Hq. In addition, note that v transforms 
as a vector under rotations within the point symmetry 
group of the lattice. 

Consider parity transformations about some point Tq, 
defined by Vra '. ro-l-r ^ tq— r. Since pseudospin is parity 
invariant {i.e. angular momentum is a pseudovector), the 
effect of a parity transformation 7'(i+j)/2 with origin at 
the midpoint between sites i and j is 

'Pii+3)/2 ■ ^ Vj., = -Vij . (7) 

Therefore, if the lattice is invariant under 7^(i+j)/2, then 

^(Vi)/2-^o^(»+i)/2 = ^0, which implies that Vy = 0, i.e. 
the hopping between sites i and j is pseudospin symmet- 
ric. In practice, the hyperkagome lattice has an intrin- 
sic handedness, and therefore is, in general, not invari- 
ant under any parity transformation. However, at low 
orders in t/Ep, many of thc hopping processes involve 
parity-invariant clusters of sites, leading to pseudospin 
symmetric hopping amplitudes. 



1. Nearest neighbor hopping 

To leading order, the nearest neighbor hopping be- 
tween sites A and C involves the creation of an electron 
at site A with pseudospin a, hopping to site C through 
the term Hac in the Hamiltonian, and projection into a 
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final pseudospin a' state: 



tca' ,Aa 



-{Ca'\HAc\Aa) 



1 



(8) 



Note that, at this level in perturbation theory, the nearest 
neighbor hopping is isotropic in pseudospin. This can be 
understood from the symmetry considerations discussed 
above - clearly, the cluster composed of the two sites A 
and C is parity invariant about its midpoint. Note that 
nearest-neighbor hopping mediated by oxygens vanishes 
to second order in perturbation theory. Thus, the leading 
correction to Eq. ([5]) occurs at third order in perturbation 
theory, and is given by 



tco 




(9) 



Even in this case, the perturbation theory processes in- 
volve parity invariant clusters, and therefore only rep- 
resent a weak pseudospin-symmetric renormalization of 
the result in Eq. ([8|) . The leading pseudospin-dependent 
terms involve next-nearest neighbor hopping, which we 
study next. 



2. Second neighbor hopping 

To third order in perturbation theory, hopping between 
sites A and D can occur through three independent paths: 

1 ^ 6 ^ D, 

2 ^ 6 ^ D, 
A^ 1 ^ 5 ^ D. 

Among these, the first path actually involves a par- 
ity symmetric cluster. However, the other two do not, 
and therefore give rise to pseudospin-dependent hopping. 
Summing the contributions from all three paths yields. 



toa' ,Aa 



where 



t 




(10) 



fll) 



The hopping amplitude from E to B can be obtained 
from Eq. pH)) through a C2 rotation at site C, which 

yields 



and for hopping from B to D, by doing an inversion 
about C, followed by a proper rotation (this operation 
is a symmetry for the clusters involved at this order in 
perturbation theory). This leads to. 



t 




(13) 



3. Third neighbor hopping 

The third neighbor hopping between sites A and E 
involves 16 different paths. The calculation is signifi- 
cantly simplified by the observation that to the lowest 
order in perturbation theory (third order), third neigh- 
bor hopping involves a reflection-symmetric cluster, and 
is therefore pseudospin symmetric. Wc find. 



tEa' .Aa 



ZEj 



(14) 



B. Parametrization of the hopping amplitudes 

Inspection of the equations P^ . and shows 
that the hopping amplitudes can be parametrized as, 

tcosO, (15) 
tsinO, (16) 
^l^.l- (17) 



V 

V 



where t represents an overall energy scale for the hopping 
Hamiltonian, 77 is the dimensionless ratio between and 
\t^\, and 9 is an angle representing the relative strength 
of direct Ir-Ir hopping and hopping mediated by oxygens. 
We expect 77 to be a small parameter. In our computa- 
tions, wc find that its precise value is not very important, 
and take r] — OA for concreteness. 

On the other hand, it is difficult to make an a priori 
estimation of the angle 6. In the case of Ir, the outer 
5d electrons are extended and there can be in princi- 
ple a sizeable orbital overlap between adjacent Ir atoms. 
Hence, we do not try to estimate the microscopic value 
of the parameter 6, which controls the relative size of 
direct Ir-Ir hopping and hopping mediated by the oxy- 
gens. Instead, in this analysis we will take 6* to be a phe- 
nomenological parameter, which will be obtained from a 
comparison with experimental data. However, as we will 
emphasize, many of the qualitative results of our analysis 
do not depend on the detailed value of 6. 



III. NON-INTERACTING LIMIT 



''pd 

iEl 



-2{U + t^) 



In the previous section, we derived a Hubbard model 
(12) for Na4lr308, see Eq. (HI). In this section, we will com- 
pute various physical observables of this model in the 
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non- interacting limit. One reason to focus on this limit 
is that experimental evidence indicates that the sys- 
tem becomes metallic under hydrostatic pressure and by 
dopingi^. Therefore, many of the properties of the metal- 
lic state may be captured by the non-interacting model, 
up to quasiparticle renormalization effects. Moreover, 
the non-interacting limit can also give us relevant infor- 
mation for the insulating state. As was argued in Ref. @, 
the insulating state is a weak Mott insulator, described 
by gapped rotors coexisting with a gapless Fermi surface 
of spinous. In this picture, the spinous would have a sim- 
ilar dispersion to the electronic dispersion obtained in the 
non-interacting model. Later on, we will discuss some of 
the effects of reintroducing interactions. 
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A. Spectrum 

The unit cell in the hyperkagome lattice consists of 12 
sites. Hence, when taking into account the pseudospin 
degree of freedom, we obtain 24 different quasiparticle 
bands in the Brillouin zone. Figure [5] shows the quasi- 
particle spectrum for the case 9 = 0, i.e. with nearest- 
neighbour hopping only. In this case, the system is pseu- 
dospin symmetric, so that pseudospin up and down bands 
are always degenerate. At the high end of the spec- 
trum, there are eight degenerate flat bands (four for each 
pseudospin), corresponding to localized states that arise 
due to the geometric frustration of the lattice. At half- 
filling, for each value of pseudospin, the chemical poten- 
tial cuts through one small electron-like pocket and two 
small hole-like pockets. Figure [3] shows the associated 
density of states. Note that near the chemical potential 
the density of states is small. This is due to the small 
size of the electron and hole pockets at the chemical po- 
tential, and leads to a strong temperature dependence of 
the specific heat in the nearest-neighbour model seen in 
Rcfs.ll audi. 

For 0^0 the hopping becomes pseudospin-depcndent. 
Figures H] and [5] show the quasiparticle spectrum and 
associated density of states for the choice of parameter 
9 = 3.31. Note that many of the degeneracies in Fig. [2] 
have been lifted, although some residual degeneracy is 
left at the high symmetry points. We have verified that 
the degeneracies at the F = (0, 0, 0) point are consistent 
with the symmetry of the magnetic point group of the 
hyperkagome lattice. Also note that, due to the further 
neighbour hopping, the flat bands seen in Fig.[2]havc now 
acquired a finite dispersion. 

At half-filling, there is a sizeable electron pocket near 
the X = (tt, 0, 0) point, and a large number of smaller 
pockets with shallow energy dispersions. These include a 
small electron pocket near (tt, 0, 0), as well as hole pockets 
near (0,0,0), and (7r/2, tt, tt). The existence of multiple 
shallow pockets leads to a strong energy dependence of 
the density of states and to a suppressed density of states 
at the chemical potential, as shown in Fig. [51 This energy 
dependence is reflected in many of the physical quantities 



FIG. 2: Spectrum along high-symmetry directions, for the 
case with nearest-neighbor hopping only 6 = 0. The momen- 
tum labels are V = (0,0,0), X = (7r,0,0), M = (7r,7r,0), 
and R = (7r,7r,7r). There are eight degenerate flat bands at 
the top of the spectrum, e = 2t. The chemical potential at 
half-fiUing is shown as a dashed line. 
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FIG. 3: Density of states for a system with nearest-neighbor 
hopping only 6 = 0. The energies are shifted by the chemical 
potential at half-filling. The delta function at high end of the 
spectrum arises from the flat bands. 

computed below. Although we only present results for 
two different values of 6*, we find that the strong energy 
dependence in the density of states seen in Fig. [3] for 
9 = 0, and in Fig. [5] for = 3.31, is generically present 
for most values of 9. This arises due to the large unit 
cell of the hyperkagome lattice, which leads to multiple 
narrow bands that straddle the chemical potential. 



B. Thermodynamics 

There arc two related experimental facts in Na4lr308 
which we will take as guides to set the parameter 9: first, 
the specific heat coefficient 7 = C/T, is very small; and 

second, the Wilson ratio W, defined hy W = ^ , is 
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FIG. 4: Spectrum for 6 — 3.31. Note that many of the de- 
generacies in Fig. [2] have been hfted by furthest neighbour 
hopping. The chemical potential at half-filling is shown as a 
dashed line. 
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FIG. 5: Density of states for 6 — 3.31. Energies measured 
relative to the chemical potential at half-filling 

much greater than one. Here, x is the low temperature 
paramagnetic susceptibility. 

Figure [5] displays the Pauli spin susceptibility the 
specific heat coefficient 7, and the Wilson ratio W, at 
low temperatures as a function of 9. For 6 = 0, the spe- 
cific heat coeflticient is very small, but the Wilson ratio 
is equal to one. This is expected, as the normalization 
of the Wilson ratio is chosen such that an isotropic non- 
interacting Fermi gas has W = 1. On the other hand, as 
6 is increased, the Wilson ratio becomes greater than one. 
However, this occurs at the price of an increased 7. At 
the value 9 = 3.31 (indicated by an arrow in Fig. |6l), we 
obtain a good compromise between a large Wilson ratio 
and a small 7. This value of 9 is in qualitative agree- 
ment with the electronic structure predicted by density 
functional theory calculationa^i. In what follows, we will 
focus on this value of 9. 

Figure [7] shows the temperature dependence of the spe- 
cific heat at 9 = 3.31. Note that C has significant cur- 
vature in the low temperature regime, which arises due 



FIG. 6: Low temperature thermodynamic quantities as a 
function of the parameter 9. Shown are the Wilson ratio 

^ ^/fef 1 the specific heat coefficient 7 = C/T, and the sus- 
ceptibility (scaled by co = Tv^k%/3fi%). The value d = 3.31 
is indicated by a vertical arrow. For concreteness, we choose 
the overall hopping amplitude to be t = 200 meV. 
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FIG. 7: Specific heat versus temperature at 6 = 3.31. Insets: 
Wilson ratio and C/T versus T for e = 3.31. C/T in the inset 
is measured in units of mj/mol K^, for t = 200 meV. 

to the strong energy dependence of the density of states 
shown in Fig. [5l The top inset in Fig. [7] shows the ra- 
tio C/T as a function of temperature, which displays a 
strong suppression at low temperatures. The bottom in- 
set shows the Wilson ratio, which increases strongly as 
temperature is reduced. This temperature dependence 
is dominated by the specific heat, as the susceptibility is 
approximately constant over this temperature range. 



C. Transport in doped samples 

Recent experiments indicate that doping Na4lr308 
renders the system metallioi^. To investigate the trans- 
port properties in such a metal, we consider the electri- 
cal and Hall conductivities in dopings ranging from 1% 
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FIG. 8: Conductivity for doping away from half- filling for 
6 = 3.3f. The conductivity is shown in arbitrary units - by 
Eq. (|18|) . (J XX is proportional to a scattering time r which is 
not known. 

to 4%, which is controlled by varying the chemical po- 
tential. We then compute transport properties using the 
Boltzmann equation results^^: 




where a = 1 ... 24 is a band index and = Vkf is the 
fermion velocity of electrons in band a. Here we have 
assumed that the scattering time t is independent of the 
wave vector k and also of the band index a. 

Fig. |5] shows electrical conductivity as a function of 
temperature, for various doping levels. Notice that for 
small doping away from half-filling, the electrical con- 
ductivity is a monotonically increasing function of tem- 
perature. This is an artifact of the assumption that the 
scattering time r used in Eq. is a constant. In other 
words, in these calculations we are only taking into ac- 
count the elastic scattering contribution to r, and ignor- 
ing the temperature dependent inelastic scattering pro- 
cesses. Thus, the temperature dependence in Fig. [8] is a 
density of states effect. More realistically, in real samples 
at high temperatures, the electrical conductivity will de- 
crease with heating. However, it is possible that at lower 
temperatures the electrical conductivity may display this 
density-of-states-induced conductivity increase. This ef- 
fect should be more pronounced for dirty samples, for 
which the elastic scattering dominates over the inelastic 
scattering over a broader range of temperatures. 

Figure [5] shows the Hall coefficient Rh as a function 
of temperature for 9 = 3.31. In contrast to the electri- 
cal conductivity, the Hall coefBcient Rh ~ — r^fj- is 
independent of scattering time, at least for wave vector- 
independent T. Note that there is a change from electron- 
like to hole-like Rh that is mediated by doping. For all 



FIG. 9: Hall coefficient for doping away from half-filling for 
e = 3.31 

dopings, we find that the Hall coefficient is a strong func- 
tion of temperature. As with other observables studied 
before, this effect is due to the presence of multiple nar- 
row bands near the chemical potential. Thus, although 
the detailed form of Rh does depend on 6 and also on 
doping, the fact that it is a strong function of tempera- 
ture is a generic feature of the model in Eq. ^ . 



IV. DISCUSSION ON THE SPIN LIQUID 
PHASE 

Although the tight-binding spectrum is derived for 
electrons in a metal, the same tight binding model would 
describe the spinon spectrum in the spin liquid phase 
up to renormalization effects coming from the interac- 
tion between spinous and the gauge field. It was shown 
that this interaction only leads to weak renormalization 
effects near the critical point between the metal and the 
spin liquid^. For example, it only leads to a Inlnl/T 
enhancement of the specific heat coefficient C/T^ On 
the other hand, the previous analysis did not include the 
spin-orbit coupling and it is expected that there would 
be an additional spin-orbit-coupling-induced renormal- 
ization of specific heat and susceptibility. That is, just 
like the metallic phase discussed earlier, the Wilson ratio 
would be in general bigger than one: the same arguments 
with the pseudospin-dependcnt hopping amplitudes ap- 
ply to the spinous in the spin liquid phase. 

There is, however, even further renormalization of the 
Wilson ratio due to the combined effect of the spin-orbit 
coupling and the quasiparticle (electrons or spinous) in- 
teractions. In the framework of the Landau Fermi liq- 
uid theory, it was shown that additional quasiparticle 
interactions that depend on both the orbital and spin 
quantum numbers can arise in spin-orbital coupled sys- 
tems. In particular, an additional Fermi liquid renor- 
malization factor for the Wilson ratio has the form: 
fi+F ifi+cN^-CG ' '^tisre Fq corresponds to the usual 



8 



spin density-spin density interaction, and Gi and G2 rep- 
resent additional spin-orbit coupling-indueed quasiparti- 
cle interactions that vanish in the absence of the spin- 
orbit couphng2^. Given that the nature of the underlying 
quasiparticle interactions among spinons is quite differ- 
ent from that between electrons, it is conceivable that 
the Wilson ratio in metal and the spin liquid can also be 
quantitatively quite different. Since the interactions be- 
tween spinons are generally believed to be much stronger, 
we may expect a bigger Wilson ratio in the spin liquid 
phase. Quantitative estimation would require the deriva- 
tion of the full Fermi liquid interaction function, which 
may be an excellent topic of future study. 
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